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ABSTRACT 

In the latest in our series of papers on XMM-Newton and ground-based optical follow-up of new 
candidate magnetic cataclysmic variables (mCVs) found in the Sloan Digital Sky Survey, we report 
classifications of three systems: SDSS J144659. 95+025330.3, SDSS J205017.84-053626.8, and SDSS 
J210131. 26+105251.5. Both the X-ray and optical fluxes of SDSS J1446+02 are modulated on a period 
of 48.7+ 0.5 min, with the X-ray modulation showing the characteristic energy dependence of photo- 
electric absorption seen in many intermediate polars (IP). A longer period modulation and radial 
velocity variation is also seen at around 4 hrs, though neither dataset is long enough to constrain 
this longer, likely orbital, period well. SDSS J2050-05 appears to be an example of the most highly 
magnetized class of mCV, a disk-less, stream-fed polar. Its 1.57 hr orbital period is well-constrained 
via optical eclipse timings; in the X-ray it shows both eclipses and an underlying strong, smooth 
modulation. In this case, broadly phase-resolved spectral fits indicate that this change in flux is the 
result of a varying normalization of the dominant component (a 41 keV MEKAL), plus the addition 
of a partial covering absorber during the lower flux interval. SDSSJ2101+10 is a more perplexing 
system to categorize: its X-ray and optical fluxes exhibit no large periodic modulations; there are 
only barely detectable changes in the velocity structure of its optical emission lines; the X-ray spectra 
require only absorption by the interstellar medium; and the temperatures of the MEKAL fits are low, 
with maximum temperature components of either 10 or 25 keV. We conclude that SDSS J2101+10 
can not be an IP, nor likely a polar, but is rather most likely a disc accretor- a low inclination SW 



Sex star. 

Subject headings: 



individual: (SDSS J144659.95+025330.3, SDSS J205017.84-053626.8, SDSS 
J210131. 26+105251. 5 ) — novae, cataclysmic variables — stars: magnetic — X- 
rays: stars 



1. INTRODUCTION 

The Sloan Digital Sky Survey rSDSS: lYo7k"et alJ2000D 

has provided a wealth of new cataclysmic variables; the 
close binary systems with active accretion fro m a late 
main sequence star to a white dwarf (reviewed in lWarnerl 
1995). Due to its sensitivity down to ~21st mag, SDSS 
is especially suited to the discovery of faint CVs with 
low accretion rates and short orbital periods. Included 
in the discoveries are dozens of systems with a notice- 
able emission line of He II, which is often a signature 
of a white dwarf with a high magnetic field. The high- 
est field systems (polars) have no accretion disk and the 
white dwarf spin is synchronized to the orbit. The inter- 
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mediate polars (IPs) usually have some outer disk, with 
magnetic curtains channelling the inner disk material to 
the magnetic poles of the white dwarf which is spinning 
faster than the orbital timescale. Polars can be identi- 
fied by circular polarization, cyclotron harmonics and/or 
Zeem an splitting (review in iWickramasinghe Sz Ferrariol 
|2000(K while IPs are found through the detection of their 
spin and orbital periods. In addition, there is a class 
of high accretion rate disk systems with possible low 
magnetic field white dwarfs l|Rodriguez- Gil et alJl200l . 
termed SW Sex stars, which can also show strong He II 
lines in their optical spectra. 

For the last few years, we have used XMM-Newton to 
identify the nature of the CVs with He II and to explore 
the nature of X-ray heating in these systems. Polars 
generally exhibit hard X-rays from the accretion shock 
and soft X-rays from the white dwarf surface heated by 
the hard X-rays. The ratio of soft to hard X-rays shows 
a dramatic change with magnetic field strength and ac- 
cretion rate. The IPs are generally hard X-ray emitters 
with large absorption effects from the accretion curtains. 
Usually the spin period is much more pronounced in the 
X-ray than the optical. The SW Sex stars usually show 
no X-ray variation and are hard X-ray emitters with low 
absorption. 

The 3 sources in this study are all about 18th 
mag in the SDSS g filter a nd were all i dentifie d 
as potential magnetic CVs in iSzkodv et alJ (|2p03a). 
SDSSJ144659.95+025330.3 (hereafter SDSS J1446+02) 
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showed no polarization and limited time-resolved 
spectra could not determine the orbital period. 
SDSSJ205017.85-053626.8 (SDSS J2050-05) has the 
strongest He II line of the 3 objects, with a peak flux 
comparable to H/3 and showed evidence of high and low 
states of mass transfer. A short span of time-resolved 
spectroscopy showed an orbital period near 2 hrs. Fol- 
lowup high speed photometry bv lWoudt et al.1 l)2004|) re- 
vealed brief (260s), deep (1.5 mag) eclipses and deter- 
mined a precise period of 1.57 hrs. The third system, 
SDSS J210131. 26+105251. 5 (SDSS J2101+10) has very 
strong Balmer lines but did not show any velocity varia- 
tion during 1.3 hr of time-resolved spectroscopy. 

The XMM- Newton data obtained on these 3 sources, 
combined with additional photometry, spectroscopy and 
spectropolarimetry has allowed us to determine that 
SDSS J1446+02 is an IP, SDSS J2050-05 is a likely polar, 
while SDSS J2101+10 escapes an easy classification, but 
appears most likely to be an SW Sex star. Our observa- 
tions and results are described below. 

2. OBSERVATIONS 
2.1. X-ray 

Table ^ summarizes our X-ray and optical observa- 
tions. SDSSJ1446+02 was observed twice by XMM- 
Newton, as the 2004 January 30 observation was badly 
affected by high charged particle background throughout, 
and the X-ray exposures were curtailed after only a few 
kiloseconds. Since the X-ray instrumentation (compris- 
ing three X-ray telescopes, b acked by the two M PS and 
one pn EPIC CCD cameras. iTurner et alJl20?)l was set 
as the prime instrument, an automatic reobservation was 
triggered. However, the 2004 optical data from X MM- 
Newton's optical monitor (OM, IStriider et al.12 001'1 were 
good throughout the scheduled ~13 ks observation. In 
the second observation on 2005 January 12, good data 
were obtained from all three EPIC spectro-imagers, but 
this time the OM only succeeded for 2 out of an ex- 
pected 4 exposures. We note that for both observations 
no useful data were av ailable from the two R eflection 
Grating Spectrographs l)den Herder et al.ll200lD arrayed 
in the optical path of the MOS detectors, as the source 
is too faint to furnish adequate signal in the dispersed 
spectra. The initial visits made to SDSS J2050-05 and 
SDSSJ2101+10 were highly successful, with only a lit- 
tle background flaring in the last 3 ks of the 12 ks ob- 
servation of SDSS J2050-05, and none affecting that of 
SDSS J2101+10. Useful data were therefore obtained 
from EPIC imaging and the OM, though once again nei- 
ther target was bright enough for the RGS. 

For the extraction of spectra and lightcurves from the 
XMM-Newton data we used the tools available in the 
Science Analysis System (SAS 9 ) version 6.5.0, with cal- 
ibration files current to 2005 December 15, and followed 
the standard protocol as given by the ESA XMM-Newton 
website and the ABC guide 10 from the US GOF. In all 
cases we produced new event list files from the Observa- 
tion Data Files to incorporate the latest calibration up- 
dates and then filtered them with the standard canned 
expressions. We defined circular source extraction re- 



9 Available from http://xmm.vilspa.csa.es/cxtcrnal/ 
xmm_sw_cal/sas.shtml 

10 http:/ /heasarc. gsfc.nasa.gov/docs/xmm/abc/abc. html 



gions (radii of 480 and 360 pixels for the MOS and pn 
camera respectively), centered on the centroid of the 
source, enclosing ~ 80% of the energy to optimize for 
signal-to-noise. Simple annular background extraction 
regions were defined on the same central MOS chip for 
those 2 cameras; for the pn we used adjacent rectangu- 
lar regions at similar detector Y locations to the target. 
We also reprocessed the OM data with omf chain extract- 
ing lightcurves using more appropriate (for our relatively 
faint targets) source aperture and background regions 
to maximize signal-to-noise. We also chose binning to 
match that of the ground-based photometry, and lastly, 
to aid comparison, converted count rates to B magni- 
tudes. 

The SAS task evselect performed the extractions of 
both X-ray spectra and lightcurves for source and back- 
ground. For SDSSJ1446+02 and SDSSJ2050-05 fur- 
ther time filtering was applied to the event lists prior 
to spectral extraction to excise intervals of high X-ray 
background; the revised good time intervals (GTIs) were 
generated by setting limits on the count rate in hard 
(> lOkeV) lightcurves for the entire detector. We also 
restricted acceptable events to singles in the pn (pat- 
tcrn=0) but up to quadruples in the MOS (pattern< 12) 
to further improve energy resolution. In contrast, in 
extracting lightcurves only the standard GTIs were in- 
voked, and we kept both singles and doubles for the pn 
(pattern < 4), to both maximize lightcurve coverage and 
signal-to-noise. For SDSSJ2101+10 we also restricted 
the energy range to below 2.5 keV, where the source 
contributions dominate. SAS tasks rmf gen, arf gen and 
backscale generated appropriate redistribution matrix 
(rmf) and ancillary response files (arf) and calculated the 
relative scaling of source to background. 

The final steps in data reduction utilized general pur- 
pose utilities in the FTOOLS 11 software suite: grppha 
grouped the spectral bins and associated various files 
ready for spectral analysis in XSPEC; lcmath created 
background subtracted lightcurves, and combined the 
two MOS lightcurves, while earth2sun applied a correc- 
tion to the time stamps for the solar system barycenter. 

2.2. Optical 

Ground based photometry was obtained for all three 
sources using the US Naval Observatory Flagstaff Sta- 
tion (NOFS) lm telescope and a 2048 x2048 CCD. An 
open filter was used which is close to a V response and 
the magnitudes were calibrated from separate nights of 
all-sky photometry with Landolt standards. Light curves 
were made by using differential photometry with respect 
to comparison stars in each field and the magnitudes 
measured using IRAF 12 routines. For SDSS J1446+02, 4 
hrs of photometry began about one hour after the XMM- 
Newton observation ended on 2004 Jan 30. No data 
were obtained for the second observation but an addi- 
tional 5 hrs of photometry from 2003 May 07 was used 
to establish the optical period. For SDSS J2050-05, 3 
hrs of photometry were obtained 3 days after the XMM- 

http:/ /heasarc. gsfc.nasa.gov/lhcasoft/ftools/ 
12 IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy 
(AURA) Inc., under cooperative agreement with the National Sci- 
ence Foundation 
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TABLE 1 

Observation Summary 



SDSS J 




UT Date 


Obs 


UT Time 


Ch.cireict.cristics a 


Comments 


SDSS J1446- 


-02 


2002 


May 11 


SO: Bok (2.3m) 






v = +0.38% 


3ks exposure, spectropolarimetry 






2003 


Apr 27 


APO 


05:43 


- 09:47 




23 x 600s spectra 






2003 


May 07 


NOFS 


06:39 


- 11:34 


~V = 17.9 


open filter photometry 






2004 


Jan 30 


XMM-Newton: 


















EPIC-pn 


04:17 


- 05:32 


0.2 cts s" 1 


only 3.0ks live time b,c 










EPIC-MOS1/2 


03:55 


- 05:18 


0.1 cts s _1 


only 4.5ks live time c 










OM 


04:03 


- 08:10 


B = 18.3 


13.9ks duration 










NOFS 


09:13 


- 13:18 


~v = 18.2 


open filter photometry 






2005 


Jan 12 


XMM-Newton: 


















EPIC-pn 


06:22 


- 08:09 


0.20 cts s _1 


5.6ks live time 










EPIC-MOS1/2 


06:00 


- 08:14 


0.10 cts s" 1 


7.7ks live time 










OM 


06:08 


- 08:16 


B = 18.2 


5.6ks duration 


SDSS J2050- 


05 


2003 


May 29 


MMT 






v = +1.06% 


1.2ks exposure, spectropolarimetry 






2003 


May 30 


MMT 






v = -0.23% 


2.4ks exposure, spectropolarimetry 






2004 


May 14 


SO: Bok (2.3m) 






v = -0.31% 


1.2ks exposure, spectropolarimetry 






2004 


Oct 18 


XMM-Newton: 


















EPIC-pn 


10:13 


- 13:17 


0.67 cts s" 1 


8.7ks live time b 










EPIC-MOS1/2 


09:50 


- 13:22 


0.21 cts s" 1 


12ks live time 










OM 


06:08 


- 08:16 


B = 18.4 


8.4ks duration 






2004 


Oct 21 


NOFS 


01:33 


- 04:34 


~V = 18.0 


open filter photometry 


SDSS J210H 


-10 


2003 


May 29 


SO 






v = +0.03% 


1.6ks exposure, spectropolarimetry 






2003 


Sep 22 


SO: Bok (2.3m) 






v = +0.10% 


1.6ks exp., i?-band imaging polarimetry 






2004 


May 14 


APO 


10:47 


- 11:04 




Iks spectrum 






2004 


May 19 


NOFS 


08:01 


- 11:35 


~V = 18.9 


open filter photometry 






2004 


May 19 


XMM-Newton: 


















EPIC-pn 


16:57 


- 18:13 


0.64 cts s" 1 


4. Iks live time b 










EPIC-MOS1/2 


16:35 


- 18:18 


0.23 cts s" 1 


6. Iks live time 










OM 


16:43 


- 18:20 


B = 19.0 


5.5ks duration 



a The open filter photometry from NOFS has an estimated equivalent V zero-point, while for spectra the flux density at ~5500A 
is used. The XMM-Newton count rates are average values for each observation for a single detector. 

b The live time of the X-ray CCD detectors refers to the sum of the good-time intervals, less any dead time. It is typically 
much less than the difference of observation start and stop times. 
c These X-ray exposures were curtailed by severe particle background flaring. 



Newton observations, while the 3.5 hrs of observations 
for SDSS J2101+10 took place 5.5 hrs prior to the start 
of XMM-Newton coverage. In all cases, the ground-based 
measurements agreed with the OM in showing the sys- 
tems were all in their normal state of accretion. 

In order to determine the orbital period for 
SDSS J1446+02, 4 hrs of time-resolved spectra were ob- 
tained on 2003 Apr 27. The Double Imaging Spec- 
trograph (DIS) was used on the 3.5m telescope of the 
Apache Point Observatory (APO). Twenty three 10 min 
blue and red spectra were obtained with a 1.5" slit, cover- 
ing the regions 4200-5100A and 6300-7200A with a res- 
olution of about 2A. IRAF routines were used to cal- 
ibrate the spectra for wavelength and flux using stan- 
dards from the night. Velocities were measured for the 
prominent emission lines using the centroid ( "e" ) routine 
in the IRAF s plot package; a double-Gaussian method 
( Shatter 1983) was tried as well. A single spectrum 
of SDSSJ2101+10 was also obtained at APO 5 nights 
before the XMM-Newton observation. The Balmer and 
He II str engths are very sim ilar to the SDSS spectrum 
shown in iSzkodv et alJ l)2003aj) . 

Circular spectropolarimetry was also performed for 
the three systems in a search for evidence of magnetic 
fields. The CCD s pectropolarimeter SPOL was used 
( Sch midt et al.llT992|) on the Steward Observatory 2.3 m 
Bok telescope and the 6.5 m MMT, as indicated in Ta- 
ble 1. All measurements on SDSS J1446+02 and SDSS 



J2101+10, and two of the three on SDSS J2050-05, 
yielded spectrum-averaged values \v\ = \V/I\ < 0.4%. 
Each result is well within 3a of zero, and thus consistent 
with being unpolarized. However, a third epoch on SDSS 
J2050— 05, obtained on 2003 May 29 in good observing 
conditions at the MMT, yielded \v\ = +1.06%. The S/N 
of these data is sufficient to reveal that the circular po- 
larization rises continuously from v ~ at ~4200 A to 
v ~ 3% beyond A = 8000 A. This and other evidence for 
polarization by cyclotron emission in SDSS J2050— 05 is 
discussed in §4.2. 

3. ANALYSIS AND RESULTS 

3.1. X-ray Spectral Fitting 

The extracted spectra were binned at >20 counts/bin, 
to facilitate the use of x 2 statistics to find the best model 
fits. For all sources the background contributed less than 
3% of the flux within the source aperture, hence simpler 
fitting of background subtracted spectra was deemed ac- 
ceptable. At the very lowest energies the calibration of 
the EPIC detectors remains uncertain; following the lat- 
est guidelines we restricted our fitting to > 0.2keV for the 
MOS and > 0.15keV for the pn. In every case, we per- 
formed joint fits to the data from all three X-ray cameras 
simultaneously, where all model parameters were fixed 
apart from the relative normalization of the pn relative 
to the MOS (which were assumed to be identical). In 
the cases where we separated the data into two phase 
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SDSS 1446 




Orbital Phase 

Fig. 1. — Radial velocities for the Ha emission line in 
SDSS J1446+02, phase folded on the best fit period, P = 3.8 ± 0.3 
hrs. Overplotted is the sinusoid fit with K = 105 ± 11 km s _1 and 
7 = —8 ± 5 km s — 1 . 



or time bins, we then had 6 datasets. Starting from the 
case where all relative normalizations, as well as other 
model parameters, were kept constant, we then allowed 
normalizations to vary (but keeping the MOS to pn ratios 
constant), and then other model parameters like temper- 
ature, F-testing to see whether the additional degrees of 
freedom were statistically warranted each time. The best 
fits we present in figures and tables are those where the 
model has the fewest free parameters (i.e. leaves the 
most degrees of freedom) ; any further freeing of param- 
eters did not then signficantly reduce the (reduced) % 2 
values. 

Furthermore, in finding the best fit models we always 
started from the simplest single emission component, 
solely absorbed by interstellar dust, with an initial value 
for Nh as predicted by dust maps for the appropriate 
line of sight through the Galaxy. Typically, we used a 
single bremsstrahlung or a version including explicit line 
emission- the XSPEC model mekal- as expected for the 
optically thin thermally-emitting plasma encountered in 
CVs. In all cases, these simplest models failed to pro- 
vide a good fit, hence we moved onto a variety of more 
complex combinations, e.g. multi-temperature or two- 
temperature thermal plasmas, additional soft blackbody 
components (as emitted by the heated polar caps in po- 
lars) and the effects of local partial covering absorption 
(i.e. due to obscuration by the accretion stream and/or 
curtain). We discuss the details of the best fits, as we 
report the results of our various analyses for each target 
in turn. 

3.2. SDSSJ1U6+02 

3.2.1. Search for periodicities 

From our APO spectroscopy run on SDSS J1446+02, 
we generated radial velocity curves for the prominent 
emission lines. A least square sine fit to the velocities 
was then used to determine the systemic velocity, the 
semi-amplitude, the orbital period and the phase (based 
on the red-to-blue crossing). While the Ha, H/3, H7 and 
He II lines were all measured, there were large deviations 
with respect to the best fit sine wave in all cases (total 
sigma of the fit of 34, 46, 75 and 73 km/s respectively), 
hence we only consider the fit to the Ha curve in any 
detail (shown in fig. [TJ. All lines gave a period solu- 
tion near 4 hrs, which is very close to the length of the 
dataset. The Ha solution shown in figure^is for a period 
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Fig. 2. — The NOFS ~V lightcurves from the two different 
observations of SDSS J1446+02 are plotted using a common time- 
scaling. We also over-plot vertical bars to mark the length of each 
spin cycle; comparison of the variations during each successive spin 
cycle shows how much more regular they are in 2003 May 07 than 
in 2004 Jan 30; possibly enhanced flickering is responsible for the 
increased scatter in the latter curve. 



of 3.8 ±0.3 hrs, K = 105 ±11 km s" 1 and 7 of -8 ±5 km 
s _1 , with red-to-blue crossing (phase 0) at 7:30 UT on 
2003 Apr 27. Because the period is so close to the length 
of the observation, it is not a robust determination, but 
useful in showing that the orbital period is not short (i.e. 
under 2 hrs). 

For SDSSJ1446±02 we also possess four optical 
lightcurves (two from the ground and two from the OM) , 
plus those from the EPIC X-ray cameras. To maxi- 
mize signal-to-noise we summed the results from the two 
MOS; however, since these are always of longer dura- 
tion than the pn, we did not combine pn and MOS at 
this stage. W e performed a Lomb-Scargle periodogram 
ijScargld Il982f) analysis of each of our lightcurves. In 
every case, a significant peak was found at a period 
of 49 min, which we identify as the spin period of the 
white dwarf. In addition, longer term variations were 
apparent in the optical lightcurves with peaks corre- 
sponding to around 4 hrs, close to the signal found in 
the time-resolved APO spectra. To determine the pe- 
riods most accurately, we used sinusoidal fitting to the 
longest 2003 May 7 NOFS lightcurve, with a sinusoid 
plus first harmonic model for the pulse, and a simple si- 
nusoid for the longer period variation. This yielded a 
period of 0.0338 ± 0.004d (48.7 ± 0.5 min) for the pulse, 
and 0.165 ± 0.015d (4.0 ± 0.4hr) for the longer period, 
which we presume to be the orbital period of the sys- 
tem. Once again we caution that the 2003 May curve 
provides only 1.3 cycles coverage of our tentative orbital 
period, and it remains poorly constrained. Nevertheless, 
we then fit the other three optical curves with sinusoids 
constrained to the pulse and orbital periods, and sub- 
tracted the latter. In fig. we show the two resulting 
NOFS lightcurves; the pulse profile is obvious in most 
cycles in the 2003 May curve, whereas in 2004 Jan there 
is far more scatter, especially away from minimum light. 
This may simply be due to varying amounts of flickering, 
perhaps indicative of small changes in accretion state. 
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TABLE 2 

X-ray Spectral Fits for SDSS J1446+02 



Model 


reduced 

x 2 

(d.o.f.) 


xlO 20 

_9 

cm 


Partial 
(xl0 22 cm~ 2 ) 


covering 

frac. 


Emission 
model parameters 


bremss 
power law 
multi-T MEKAL b 

2 T MEKAL 


1.5 (139) 
1.11 (136) 
1.14 (137) 

1.10(134) 


6.4 a 
1 4+0-8 

1.6 ±0.7 

4±f 


0.8 ±0.2 
1.3 ±0.2 

n q+3-9 
u.y_ 2 


41+0- 04 
0.64 ± 0.02 

n 47+0-05 


kT = 200 (pegged) 

T = 1.2 ±0.2 

fcT max = 60 keV (fixed c ) 

a = 1.0 (fixed) 

fcTi = 130± 2 eV 

kT 2 = 80t% keV 




Joint fit to phase-selected 


"maximum" and "minimum" flux intervals: 




2 T MEKAL 


1.08 (104) 


13+ 2 


none (max) 
4.4±J;| (min) 


none (max) 
0.61 ±0.03(min) 


fcTi = 81±J eV 
fcT 2 = 58t 22 keV 



a This fit did not converge well, hence no error estimates are available. 

b We used the XSPEC model CEMEKL, in which a fine grid of mekal models are coadded, with the emission measures following 
a power-law in temperature with index a and up to feT max (i.e. normalizations scale as (T/Tmax)°). 
c fcT max fixed to plausible value, free fit does not converge. 



Lastly, we phase folded (and binned) the light curves on 
the pulse period to examine its profile in greater detail. 
The binned results are shown (along with the X-ray) in 
Fig. The profile is clearly far from sinusoidal, having 
a relatively broad and flat-topped peak, and narrower 
V-shapcd minimum. 

Before phase-folding the X-ray data, we created five 
different energy selected lightcurves, in order to investi- 
gate any dependence of the pulse profile on energy. The 
final binned lightcurves shown in Fig. [21 were constructed 
by summing the folded and binned results for all three 
X-ray cameras to maximize signal-to-noise. Perhaps not 
too surprising, the profile of the X-ray pulse is very dif- 
ferent from the optical, but it is noteworthy that the 
X-ray and optical are neither in phase nor anti-phased. 
Instead, the peak of the X-ray roughly leads that of the 
optical by ~0.3 in phase. The X-ray pulse is also energy 
dependent, with the largest amplitude (and symmetrical 
peaked pulse) in the two lowest energy bins, the ampli- 
tude then decreases with energy, until above 5 keV there 
is no significant variation. 

3.2.2. Spectral variations with phase 

The variation in pulse amplitude with energy is ex- 
actly that which we would expect if obscuration by a 
local absorber is responsible. In fitting the X-ray spec- 
tra, we first found a fit to the entire dataset, but then 
applied this model separately to spectra phase-selected 
from the maximum (</> = 0.05 — 0.55) and minimum 
(</> = 0.55 — 1.05) flux intervals of the pulse profile. The 
details of the fits are given in Table [3 The model that 
has the best fit consists of a two-temperature thermal 
plasma (mekal), absorbed both by a small Galactic col- 
umn, but also by a much larger variable partial covering 
absorber. The fit to the combined spectrum achieved a 
barely acceptable reduced x 2 = 1.1 for 134 degrees of 
freedom (d.o.f.), though it does better than any other 
model. This model has a well-constrained cool plasma 
of kT = 130l^eV contributing unresolved soft emission 
line structure plus a very hot component whose temper- 



ature pegged at 80 keV, the upper limit for the mekal 
model. Furthermore, the values of the Galactic and local 
columns were very poorly constrained. This is probably 
an effect of trying to fit a highly phase- variable spectrum 
with a single model (even this fairly complex one). In- 
deed, after separating the data into "max" and "min" 
subsets, a joint fit to these spectra was a significant 
improvement (see fig. 0J. There was not a significant 
change in the reduced x 2 , but we were now able to find 
well-constrained values for the absorbing columns and a 
physically plausible temperature for the hottest mekal 
component of kT = 58l 22 ,keV. Once again the fit re- 
quired a very cool mekal to account for line structure in 
the 0.3-0.4 keV range, plus O VII emission at 0.56 keV; 
indeed this parameter now pegged at its minimum value 
of 81 eV, indicating that some problems still remain (but 
the need for a very cool plasma is less problematic than 
an unphysically hot one, given the physical conditions in 
the expected complex multi-temperature accretion col- 
umn). In this fit, the partial covering absorber parame- 
ters for the "max" tended to zero, hence we fixed them 
as such; the "min" interval as expected requires a thick 
column around 4 x 10 22 cm -2 , with 61% covering frac- 
tion. In our joint fit we found no statistical support for 
varying any of the normalization parameters; the change 
between "max" and "min" is entirely accounted for by 
the introduction of the thick partial covering absorber. 
Therefore, the fully unabsorbed flux (0.01-10 keV) re- 
mains constant at 3 x 10 -11 ergcm~ 2 s _1 , whereas the 
affect of the absorption leads to clear energy-dependent 
changes in the observed fluxes: a drop of 60% (4.6 — 1.8 x 
10~ 13 erg cm -2 s^ 1 ) in the softest 0.2-2.0 keV band; a de- 
crease of only 30% (7.6 — 5.6 x 10~ 13 erg cm -2 s _1 ) in the 
2.0-5.0 keV; and no significant change in the flux above 
5 keV. 
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2005 Jan 12 
5.0-10.0 keV 




0.5 1 1.5 

Pulse phase (P=48.5 min) 



Fig. 3. — White Dwarf spin phase folded, and binned light 
curves for SDSS J1446+02. Upper: the energy-selected EPIC X- 
ray light curves; note the marked change in amplitude, increasing 
with decreasing energy. Lower: the optical light curves from the 
two XMM-Newton/OM ( in B) and two NOFS observations (in 
an approximate V band); in the NOFS light curves the change in 
amplitude, due to the variation in scatter about the mean curve 
is made very apparent. Note: (i) the simultaneous X-ray and OM 
2005 Jan 12 lightcurves share a common ephemeris, arbitarily cho- 
sen to place optical minimum at phase=0.0, the phasing of each 
of the other optical lightcurves has simply been chosen to align 
their optical minima; (ii) for clarity each of these lightcurves has 
been offset vertically, with 2 normalized X-ray counts and 0.4 mags 
respectively. 

3.3. SDSS 32050-05 

The eclipses see n in the optical photometry of 
iWoudt et alJ (|2004f) are also very apparent in the X- 
ray lightcurves, and in our B-band OM results. Their 
ephemeris is sufficiently precise that there is merely a 
0.01 cycle uncertainty at the time of our XMM-Newton 



observations. Adding our own eclipse center measure- 
ments, we are able to even further refine the ephemeris 
to: 

HJD min (TT) = 2453296.29816(6) + 0.06542463(1) *E d 

where the numbers in parenthesis indicate the la un- 
certainty in the last digit. Also, note that the time 
is given in Terrestrial time (UT=TT+64.184 s at the 
present epoch). 

Besides the eclipses, the X-ray lightcurves exhibit large 
amplitude but apparently irregular flaring behaviour. 
We divided the X-ray data into two energy bins, above 
and below 1.6 keV, roughly encompassing the same count 
rates, then phase-folded the lightcurves and we plot them 
in fig. [5] where we also plot each cycle of data with a dif- 
ferent symbol (and color in the electronic edition). Al- 
though this energy division is not physically motivated, 
it does roughly separate the energies severely affected by 
photo-electric absorption from those little affected (as 
can been seen in fig. EJ). We note, however, that the 
soft band does include contributions from physically dis- 
tinct X-ray emission regions- the heated white dwarf sur- 
face, and the accretion column- which should be borne 
in mind in any interpretations of the lightcurves. We see 
that the flaring does not repeat in each orbit, and that 
it is far more significant below 1.6 keV. The over-plotted 
binned data (stepped line) bring out the average orbital 
modulations more clearly. At the higher energies, there 
appears to be a broad peak, cut unevenly by the eclipse, 
whereas at lower energies, the rise in flux to peak does 
not occur until phase 0. As we did for SDSS J1446+02, 
we undertook X-ray spectral fitting to both the complete 
dataset, and two phase-selected intervals, here dubbed 
"peak" and "trough" , excluding the eclipse phase: these 
regions are also indicated in fig. [3J 

For the complete dataset acceptable fits (see Table |2J 
were found for two component models with Galactic 
absorption consistent with the maximum line-of-sight 
column, with a optically thin thermal component (hot 
bremsstrahlung or mekal equally) plus a soft black- 
body. However, the temperature of the hotter compo- 
nent at ~80 keV is higher than typically found in CVs. 
We then investigated the phase-selected spectra seeking 
a single model requiring the fewest varying parameters. 
The final model consists of a two-temperature thermal 
component plus a soft blackbody, with partial covering 
absorption. In the light of this success, we also fit the 
complete dataset with the same model; in all cases we 
found physically plausible temperatures and absorption 
columns. The dominant contribution to the observed 
absorbed flux (98% in 0.01-10 keV range) is from a hot 
MEKAL with kT« 40 keV, but F— tests confirm the im- 
portance of the soft blackbody to provide the softest X- 
ray flux (at 93% confidence), and the cool mekal to ac- 
count for significant unresolved line emission at ~lkeV 
(98%). The change in the spectrum from peak to trough 
is in part due to increased partial covering absorption: 
the peak spectrum requires no absorption in excess of 
the Galactic column, while for the trough we find a 35% 
covering column with Nh = 1.3 x 10 22 cm~ 2 ; and also 
to a decrease by 25% in the normalization of the hot 
mekal; the blackbody and cool mekal normalizations 
remain unchanged. The fully unabsorbed 0.01-10 keV 
flux amounts to 5 x 10~ 12 ergcm~ 2 s _1 , with 40 percent 
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Fig. 4. — Upper panel: XMM-Newton/EPlC-pn spectra for SDSS J1446+02 from the "max" (green or black points) and "min" (red or 
dark grey) spin phase intervals; the MOS data, used in the full joint fits, have been omitted for clarity. The best fit two temperature MEKAL 
model is over-plotted for each (solid and dashed respectively); we also plot the two component contributions. Lower panel: the residuals 
to the fit of the total model. 



TABLE 3 

X-ray Spectral Fits for SDSS J2050-05 



Model 


reduced 

x 2 

(d.o.f.) 


N H 
xlO 20 

_9 

cm 


Partial covering 
Njj frac. 
(xl0 22 cm" 2 ) 


fcT BB 
(eV) 


fcT(s) 
(keV) 


bremss 
MEKAL 
2 T MEKAL 

BB + MEKAL 

BB + bremss. 

2 T MEKAL + BB 


0.99 (418) 
1.01 (418) 
1.00(416) 

0.99 (416) 
0.98 (416) 
0.99 (412) 


4.3 ±0.3 
4.8 ± 0.3 

7.4 ±0.3 

0,1 -OA 
5 fi+°' 4 

°- D -o.s 

c 1+0.5 
°- 1 -0.7 


1 4+0-8 14 +0.45 
1- -0.6 u,± -0.05 


20+J 

22+ 1 
z -10 

32±J 


80l« 

(8.08+q 1 ) x 10~ 2 

8oi; 7 

801° 
100+3 2 
1 3+ 1 - 3 
59+ 21 






Joint fit to ph 


use-selected "peak" and "trough" intervals: 






2 T MEKAL + BB 


0.88 (417) 


4 9+H 


none (pk) none (pk) 
1.3lg;|(tr) 0.35 ± +0.03(tr) 




1 0+ 03 
1 - u -0.2 

41+5 



contributed by the soft blackbody. In fig. H3 we show the 
fits to the phase-selected spectra, together with a break- 
down of the three components; for clarity we only show 
data from the pn, though we used pn and the two MOS 
spectra when fitting. 



3.4. SDSSJ2101+10 

With no previously known periods, either spin or or- 
bital, our first step was to examine the lightcurves (see 
fig. 01 for periodicities. Unfortunately, the XMM-Newton 
data span only ^1.5 hr, and even the NOFS lightcurve is 
only 3.5 hrs long, limiting our sensitivity to any modula- 
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Fig. 5. — Energy-selected XMM-Newton/EPlC-MOS (combined) X-ray lightcurves of SDSS J2050-05, phase-folded on our revised 
ephemeris: (top panel) a hard band with E > 1.6keV; (bottom panel) soft band with E < 1.6 keV. Each of the three cycles covered is shown 
with a different symbol and grey shade (color in electronic edition) in order to show how much of the variability is apparently random 
flickering, and is more pronounced in the soft energy band. The average orbital modulation is brought out by binning the data; this curve 
is over-plotted in solid black. The vertical grey bars and labelling indicate the three phase intervals into which the data were subdivided 
to study variations in the spectra. 



tions with periods much in excess of 2 hours. Moreover, 
the variability in both X-ray and optical bands appears 
extremely complex. The longest, NOFS optical, curve 
appears to show two distinct humps at around 0.87 and 
0.94 (truncated HJD(TT)), but of course even here only 
two cycles are pr esent. Running a phase dispersion mini- 
mization (PDM, Stcllmgwcrl ll978fl period- folding search 
on this light curve, we find three signals, but the low- 



est frequency one is close to the inverse of the data time 
span, and is therefore unreliable. This leaves minima at 
around 13.5 d _1 (107 min) and 26. 4d -1 (55 min), which 
are close to being harmonically related to each other. In 
figure 00 we show the PDM for the NOFS lightcurve and 
the lightcurve folded on the two candidate periodicities. 
The fold on the longer period confirms the repeatability 
of the aforementioned humps: only at phases (arbitrary) 
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Energy (keV) 

Fig. 6. — Upper panel:XMM-NewtonfEPIC-pn spectra for SDSS J2050-05 from the "peak" (black or green points) and "trough" (red or 
dark grey) orbital phase intervals; the MOS data, used in the full joint fits, have been omitted for clarity. The best fit two temperature 
MEKAL plus blackbody model is over-plotted for each (solid and dashed respectively); we also plot the various component contributions. 
The hotter emitting gas provides the dominant contribution, and the change in overall flux between the intervals is due to a combination 
of its changing normalization and the extra absorption from a dense partial covering component. Lower panel: the residuals to the fit of 
the total model. 



0.25-0.5 do the two cycles agree; but the fold is also some- 
what double-humped, with possibly a second interval of 
higher flux at around 0.9, although here there is a large 
amount of scatter. Folding on 54 min the points line 
up along an approximately sinusoidal curve, but again 
there is significant scatter about this mean at all phases. 
Furthermore, for this shorter period we can usefully in- 
terrogate the OM and MOS light curves: their PDMs are 
also plotted in fig. El no strong signals appear around 54 
min. 

Returning to fig-H in order to consider the NOFS pe- 
riods further we have over-plotted vertical bars to mark 
phase 0.0 and 0.5 of a 107 min period, arranging the 
time axes for the XMM-Newton and NOFS observations 
to yield phase agreement (only 4 cycles elapse between 
the two sets of observations) . We see that the repeatable 
NOFS humps agree in phase and approximate form with 
a hump in the OM lightcurve. However, at other phases 
agreement is lost, in particular another peak precedes the 
major hump in the OM, possibly correlated with a similar 
X-ray peak in fact. Ignoring the high frequency flicker- 
ing in the X-ray lightcurves, there is possibly a slower 
variation on the *~100 min timescale, though of course 
we have no way to test whether this behaviour actually 
repeats. 

In the absence of any decisive periods, we simply fit 



models to the complete X-ray spectrum. These are de- 
tailed in Table H A range of emission components from 
a single fairly cool mekal (kT = 9keV), to a multi- 
temperature version with fcr max = 25keV provide ade- 
quate fits, all with a single absorbing component with 
Nh ~ 8 x 10 20 cm -2 consistent with the Galactic line-of- 
sight maximum. An example fit to the spectra is shown 
in fig. |5| the resulting unabsorbed 0.01-10 keV flux is 
1.7 x 10 -12 erg cm -2 s -1 . 

4. CLASSIFYING THE THREE SOURCES 

4.1. Nature of SDSS J 1446+02 

All the observational evidence indicates that 
SDSS J1446+02 is a short period IP. A classic sig- 
nature of IPs is the presence of multiple periodicities 
in their lightcurves. For SDSS J 1446+02 we find that 
both its optical and X-ray fluxes are strongly modulated 
at a short 49 min period, and there is evidence for 
an additional optical modulation at 4 hrs. This latter 
periodicity also appears in the Ha radial velocity curve. 
However, more extensive photometric and spectroscopic 
runs are needed to confirm these, and hence secure the 
determination of the orbital period. The identification 
of the 49 min modulation with the spin period of an 
asynchronously rotating white dwarf is much more se- 
cure: the characteristic energy dependence of the X-ray 
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TABLE 4 

X-ray Spectral Fits for SDSS J2101+10 



Model 


reduced 




kTx or 


kT 2 (keV) 




X 2 (d.o.f.) 


Xl0 20 cm- 2 


fcT max (keV) 


or a 


MEKAL 


0.82 (95) 


7.1 (fixed) a 


9± 1 




MEKAL 


0.80 (94) 


8 4+ 1 ' 4 


8± 2 




2 T MEKAL 


0.77 (92) 


8.6+1- 2 


™±l:l 


1.0 ±0.2 


multi-T MEKAL 


0.75 (93) 


s.6± ;i 


25±i 


1 2 +0 ' 4 



a The absorbing column is fixed to that estimated from dust maps, given by the FTOOL nH 




1.2 

HJD(TT)-2453144 



Fig. 7.— Lightcurves for SDSS J2101+10 from NOFS (optical) 
and XMM-Newton (optical and X-ray). The lower three panels 
share a common time axis, while the earlier NOFS observation 
differs and is marked on the top axis. The vertical bars mark 
the times of (arbitrary) phase 0.0 and 0.5 for a 107 min period 
(detected in the NOFS lightcurve). 



flux modulation, and our simple phase-selected spectral 
fitting confirm that the variation in flux is dominated by 
changes in the local absorbing column. This is exactly 
as seen in many other IPs, as a result of the extended 
accretion curtain intersecting our linc-of-site to the 
X-ray emitting accretion column. 

The underlying X-ray emission can be fit with that 
emitted by a two-temperature optically thin plasma, i.e. 
bremsstrahlung. No doubt this simple model only ap- 
proximates what must be a complex multi-temperature 
shock region in the accretion column, however the hotter 
component with kT sa 60keV indicates the upper limit 



to the temperatures therein; again this value is typical 
of IPs. 

Assuming that P or b = 4hr, we find a ratio of 
fspin/^orb = 0.2, placing SDSS J1446+02 well-within the 
group of so-called conventional IPs, as defined by 0.25 > 
fspin/^orb > 0.01 and P orb > 4hr <|Norton et al J 1200^1 . 
In all respects the observed properties of SDSS J 1446+02 
are consistent with a garden- variety IP. 

4.2. Nature of SDSS 32050-05 

The observational evidence favors a polar classifica- 
tion for SDSS J2050-05. One spectropolarimetric ob- 
servation found significant circular polarization that in- 
creases to v ~ 3% for A > 8000 A. This result sug- 
gests that the strongest cyclotron harmonics may lie in 
the near-IR, and thus that the magnetic field is rela- 
tively low, B < 30 MG, or that the specific accretion 
rate in the cyclotron-emitting portions of the funnel(s) is 
relatively low, m < 1 g cm~ 2 , or both. No cyclotron 
harmonics are evident in either the total flux or cir- 
cular polarization spectra, consistent with a high tem- 
perature shock. Therefore, SDSS J2050— 05 is proba- 
bly not a low-accret i on rate polar ( Schwoue et al]l2002l 



Schwop 
2005), 



ISzkodv et alJl20n5fl ISchmidt et al .11200,51). where mass- 
transfer appears to occur not via Roche-lobe overflow, 
but rather by efficient magnetic capture of the secondary 
star's stellar wind. 

In support of the above interpretation, the X-ray spec- 
tra of SDSS J2050— 05 are well-fit by a two-temperature 
thermal plasma model (mekal), but in addition the data 
require a soft blackbody component. In this case, the two 
mekal components have temperatures of 1 and 41 keV 
respectively, once again likely indicating the limits of the 
range of temperatures present in the optically thin emit- 
ting plasma, both quite typical for polars. The black- 
body has a temperature of about 30 eV, again within 
the observed range for polars. In the "peak" region, 
this cool blackbody component, arising from the heated 
surface of the white dwarf, contributes about 40 per- 
cent of the unabsorbed 0.01-10 keV flux. In terms of 
the energy balance, assuming an X-ray albedo for the 
white dwarf surface of ax = 0.3 and that the cyclotron 
contribution to the energy losses is negligible, one finds 
L BB /L hr » 7t/bb(1 - a x )d 2 /2irf hl .(l + a x )d 2 « 0.2, 
which given uncertainties in estimating the unabsorbed 
soft flux is comparable to recent measures of the energy 
balance in other actively accreting polars. 

As an eclipsing system, the orbital period is very 
well-established. In the X-ray lightcurves we again see 
eclipses, but the remainder of the variability is complex. 
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Fig. 8. — Phase Dispersion Minimization periodograms for 
the MOS X-ray and NOFS and OM optical lightcurves of 
SDSS J2101+10. The successive curves have been offset vertically 
by 0.5 for clarity, and have been truncated at the low frequency 
end at the inverse of the respective lightcurve time span. The inset 
shows the NOFS lightcurve folded on its candidate periods, 54 and 
107 min. The ordinate is in relative V magnitudes, with a 0.5 mag 
interval between large tick marks. 



It is dominated by aperiodic flickering behavior, which is 
strongest in the softer X-ray flux (below 1.6 keV). In an 
attempt to discover any underlying modulation on the or- 
bital period, we created phase-binned light curves averag- 
ing over the three cycles of data. What remains appears 
quite typical of a polar, in which one accretion pole is 
always visible (i.e. no white dwarf self-eclipse). This pro- 
duces emission modulated on the orbitally-synchronizcd 
white dwarf spin period, since the projected area of the 
emitting region changes with phase. The fact that only 
the normalization of the hotter thermal emission com- 
ponent appears to vary between "peak" and "trough" 
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Fig. 9. — XMM-Newton/EPIC spectra (points) and best fit 
model (line) for SDSS J2101+10: 6/acfc-pn; dark grey-MOSl; light 
grey-M.OS2. This model consists of a multi-temperature optically 
thin plasma model (cemekl), absorbed by N H = 7.5 X 10 19 cm -2 . 
In the lower panel the residuals to the fit are plotted. 



spectra is probably simply owing to poor statistics for the 
much smaller contributions from the cooler mekal and 
soft blackbody components. The shape of the average, 
non-eclipse harder X-ray (E > 1.6keV) lightcurve is in- 
deed approximately sinusoidal, whereas in the soft band 
the lightcurve does not rise until after the eclipse, i.e. its 
flux remains low in the phase range 0« 0.7 — 0.95. The 
energy-dependence of this feature suggests that photo- 
electric absorption may be the cause. Indeed, the spec- 
trum for the <f> k 0.4 — 0.95 "trough" phase interval 
does require a fairly high column, which partially covers 
the X-ray emission. Such pre-eclipse energy-dependent 
dipping structure is seen in many other eclipsing po- 
lars e .g. EP Dra, HU Aq r and UZ For lIRamsav et alJ 
l20?)l ISchwone et al]l200tt ISirk fc Howelllll99St) . In the 
high-quality soft X-ray lightcurves of HU Aqr, there is 
one broad dip centered on phase 0.7 and a narrower fea- 
ture at around 0.9. These are interpreted as arising from 
the accretion column obscuring the nearby heated pho- 
tosphere of the white dwarf, while the more distant ac- 
cretion stream accounts for the sharper dip, indeed the 
width and azimuth of this dip can be related to the ge- 
ometry of the region where the stream threads onto the 
white dwarf's magnetic field. Unfortunately, given the 
quality of our own data (both lacking the high signal-to- 
noise and afflicted by the flaring) we cannot identify the 
dipping components, and hence probe the accretion flow 
in greater detail; and of course from our spectral fitting 
we know that the soft flux arises both from the white 
dwarf and various parts of the accretion column, which 
would further complicate a more detailed analysis. 

4.3. Nature of SDSS 32101+10 

Although originally identified as a candidate mCV, 
given the strength of its high-excitation optical emission 
lines, the follow-up observations of SDSS J2101+10 we 
have presented do not seem to confirm this classification. 
We find no strong modulation of its X-ray flux, quite un- 
like any polar, for instance SDSS J2050-05. The X-ray 
spectrum also does not require any absorption in excess 
of the estimated Galactic column, effectively ruling out 
a typical IP. The X-ray emission is describable by a va- 
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riety of thermal plasma (mekal) models, though with 
a noticeably lower maximum temperature (kT ps 10 or 
at most 25 keV) than is found in most mCVs or indeed 
in our fits to SDSSJ1446+02 or SDSS J2050-05. This 
cooler plasma is in fact far more typical of disk-accreting 
systems. Only the SW Sex stars show He II with strength 
comparable to the magnetics; both the lack of absorption 
in the line profiles and their barely detectable radial ve- 
locities could be accounted for by a low inclination to our 
line-of-sight. 

5. CONCLUSIONS 

We have reported on observations allowing classifica- 
tion of three new SDSS CVs showing prominent He II: 
SDSSJ1446+02, SDSSJ2050-05 and SDSS J2101+10. 
Each represents a different class of CV, known to exhibit 
such high excitation emission lines. 

SDSS J1446+02 is a clear example of an IP. We find 
two periodicities at about 4 hr and 49 min, arising re- 
spectively from the orbit and asynchronously spinning 
white dwarf. The X-ray spectrum comprises a multi- 
temperature thermal plasma, with a maximum of 60 keV, 
modulated by phase-dependent local absorption. This 
likely originates from obscuration by the accretion cur- 
tains formed by the accreting material as it flows between 
the disrupted disc and the magnetic poles of the white 
dwarf. 

SDSS J2050-05 represents the most highly magnetized 
of the three, a fully spin-orbit synchronized polar, with 
an orbital period of 1.57 hr. In this case no disk forms, 
and all accretion is channeled along the field lines to 
a single dominant pole. Owing to its high inclination, 
there is a phase interval when the accretion stream ob- 
scures the X-ray emitting region, but most of the mod- 
ulation on the spin/orbit appears consistent with sim- 
ply its changing aspect. The X-ray emission has two 
components: the dominant one is characteristic of the 
post-shock bremsstrahlung-cooled optically thin plasma 
of the accretion column, here with temperatures ranging 
from 1 keV to 40 keV; this also heats the surface of the 
white dwarf leading to additional soft blackbody (30 cV) 
emission. 

In SDSS J2101+10, we suggest that the combination 
of high accretion rate and lower magnetic field strength 
allows disk-mediated accretion, but also the formation of 
He II emission lines in the optical spectra, i.e. if it were 
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at higher inclination it would appear as a classic SW Sex 
star. Its X-ray spectra indicate a rather cool thermal 
plasma (~10 keV) more typical of disk systems, and an 
unobscured line of sight, inconsistent with an IP. From 
only weak modulation (unlike a polar) of its X-ray and 
optical lightcurves we suggest an orbital period in the 
100 min range. 
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